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ABSTRACT 

The degree of influence that environment and mass have on the steUar populations 
of early-type galaxies is uncertain. In this paper we present the results of a spectro- 
scopic analysis of the stellar populations of early-type galaxies aimed at addressing 
this question. The sample of galaxies is drawn from four clusters, with {z) = 0.04, and 
their surrounding structure extending to '~ 10i?vir- We find that the distributions of 
the absorption-line strengths and the stellar population parameters age, metallicity 
and a-element abundance ratio do not differ significantly between the clusters and 
their outskirts, but the tight correlations found between these quantities and velocity 
dispersion within the clusters are weaker in their outskirts. All three stellar population 
parameters of cluster galaxies are positively correlated with velocity dispersion. Galax- 
ies in clusters form a homogeneous class of objects that have similar distributions of 
line-strengths and stellar population parameters, and follow similar scaling relations 
regardless of cluster richness or morphology. We estimate the intrinsic scatter of the 
Gaussian distribution of metallicities to be 0.3 dex, while that of the a-element abun- 
dance ratio is 0.07 dex. The e-folding time of the exponential distribution of galaxy 
ages is estimated to be 900 Myr. The intrinsic scatters of the metallicity and a-element 
abundance ratio distributions can almost entirely be accounted for by the correlations 
with velocity dispersion and the intrinsic scatter about these relations. This implies 
that a galaxies mass plays the major role in determining its stellar population. 

Key v^rords: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: stellar content - galaxies: formation. 



1 INTRODUCTION 



Tlie classical model of galaxy formation, the 
lithic collapse model jEggen et al. 1963: Larsonl 
19751: iT inslev fc GunnI Il976l : lArimoto fc YoshiT 



1974 



1987 



Bressan et al. 19941). proposes that galaxies form in a single 



massive collapse at high redshift and that subsequent 
evolution is purely passive. This model has been overtaken 
by the current dominant mo del of galaxy formation, the 
hiera r chical merging model ('T oomrtJ Il977l : ISearle fc Zimil 
1 19781 : IWhite fc Ree3 lT978l. which proposes that galaxies 
are built up through mergers: small objects form first and 
undergo a series of mergers that build up more massive 
objects. A third scenario, the revised monolithic collapse 
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model (e.g. iMerlin fc Chiosil [20061 ') . proposes that galaxies 
form in a number of rapid mergers at high redshift before 
evolving passively. 

Galaxy clusters originate from the most extreme den- 
sity fluctuations, where galaxy formation and evolution is 
expected to proceed at an accelerated rate. Moreover, the 
stellar populations of galaxies moving into a cluster envi- 
ronment are expected to be modified via interactions with 
other galaxies and the dense intra-cluster medium. These 
interactions can be roughly divide d into two bro ad classes: 
local process es including mer gers ('To omre|[l977l ) and tidal 
interactions l|Mastropietro et al. 20 05); and global proc esses 
including ram pressure stripping (iGunn fc Gottlll972l ). m- 
teracti ons with the clust er tidal field ( Bekkilll999l). harass- 
ment I Moore et al.|[l999l ). and strangulation I Larson et al.l 
Il98d ). Local processes are more efficient in galaxy groups. 
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where relative velocities are lower, while global processes are 
more efficient in clusters where the frequency of interaction 
is higher. 

Correlations are therefore e xpected b e tween 
galaxy observables and en vironment ()Kauffmannl 1 19961 : 
iKauffmann fc Chariot _ 19981) . and have b een found for 
galaxy morpholog y (IDayis fc Gelleil Il976l : iDressleJ Il980l : 



IPostman fc GeUerl Il984l : 
llBlanton et alj 12005^ 
I2OO5I: [Hashimoto fc Oemlei 



Baloeh et alT 1 19981) . colour 



Ser sic index ( Blanton et al.l 
1999|}^ star-formation rate 



dLewis et al.l I2OO2I : IComez et al.1 I2OO3I: ^elli fc Gavazzil 
l2006l ). and spectral type (|Norberg et al.ll2002l ). However, 
in a comprehensive study of the dependence of galaxy ob- 
servab les on environment within the SDSS, iBlanton et al.l 
iiooi) suggest that the structural properties of galaxies 
are less dependent on environment than their masses and 
star-formation histories. 

To add to this debate. Fundamental Plane studies have 
found no differences between field galaxies more massive 
than 2 x 10^^ and their counterpar t s in clusters at the 
same redshift (iTreu et al.1 Il999l . I2001I : Ivan Dokkum et al.1 
l200lh . while less massive galaxi es were found to be younger 
in the field (iTreu et al.H 2002: van Dokkum fc Stanfordll2003l : 
Ivan der Wei et al.1 12004 . .2005. : ITreu et al.1 l2005al lbl). This 
might imply that star formation in field galaxies occurs first 
in the most massive galaxies then progressively in less mas- 
sive galaxies, and that mass rather than environment gov- 
erns the overall growth. 

The Lick system of abs orption-line indices 
llBurstein et al.l[l984 : Fabc r et al.ll98 5: Burstoin ct al. 1986,; 
ICorgas et aL[ll993l : [Worthev et a l. 1994; Trage r et al.ill99i) 
and associated models (e.g. Iworthev. .1994 : Thomas et al.l 
I2OO3I . I2OO4I ). from which the stellar population parameters 
(SPPs) age, metallicity ([Z/H]), and a-element abundance 
ratio ([a/Fe]) can be estimated, provide an excellent method 
with which to study the role of mass and environment in 
determining galaxy properties. 

The line-strengths of early-type galaxies are found to 
be correlated with velocity dispersion (ct). Lick indices that 
are more sensiti ve to [Z/H] effects, e.g. Mg&, are posi- 
tively correlated llBurstein et al.lll988l: iBender et al.l[l99^: 
Ziegler fc Bender] Il997l : ICoUess et al.1 Il999l : iKuntschiierl 
2000h while those indices more sensitive to age effects, 



e.g. H/?, a r e neg a tively c orre lated (Fisher ct al. "l9 95l. Il996l; 



J0rgensenl 119971. 11999a 



20031 : iBernardT et al lEoo 



Jcimtschner .200Q : 
3). Although both ( 



Caldwell et al.l 
elements are sen- 
sitive to [Z/H], the slope of the correlation between Mg and 
velo city dispersion is found to be much steeper than that of 



Fe (IWorthev et all 



s lound to be mucn steeper tnan tnat 01 
I992I; iFisher et alii 19951; ICreggid [l997l; 



Kuntschneij |200GI : iTerlevich fc Forbes! 



J0rgensenl 1999al : 
20021 ). 

These index-a correlations suggest that the SPPs 
should also be correlated with velocity dispersion. Given 



the difference in the slopes of the M; 
tions, a correlation between [a/Fe] anc 
was expected and has been confirmed ( 


y-a and Fe-a rela- 
velocity dispersion 
Trager et al.l 2000a: 


Proctor fc Sansom 2002: Thomas et al.l 


2002:lMehlert et al.1 


.2003; Thomas et al. 2005). If the a-element enhancement is 



due to the timescale of star formation, and velocity disper- 
sion is a proxy for mass, then this correlation indicates that 
more massive galaxies form their stars on shorter timescales 
than less massive galaxies. 



The [Z/H] of a galaxy is also found to correlate w ith ve- 
locity dispersion (|Greggid[T997l : iThomas et al.ll2005l ). with 
massive galaxies being more metal-rich. Such a relation 
is a natural conseque nce of galactic wind models (e.g. 
lArimoto fc Yoshiilll987l ). which show that the larger gravi- 
tational potential of massive galaxies allows them to better 
retain their heavy elements. 

Whether a correlation exists between age and veloc- 
ity dispersion is still uncertain. Early studies found no sig- 
nificant correlation (I Trager et aL I l2000al : iKuntschner" et al.l 
I2OOII : ITerlevich fc Forbe^ l2002h . but recent studies have 
detected a weak but significant c orrelation hav i ng lar ge 
scatter (iProctor fc SansomI |2002| : iProctor et al] I2004al 0: 
iThomas et al.l |2005| ). with more massive galaxies being 
older. This trend, which would indicate that more massive 
galaxies formed their stellar content earlier, is in agreement 
with the concept of down-sizing (Cowie et al. 1996) and 
the latest semi-analy tic models of galaxy formation (e.g., 
iDe Lucia et~alll2006l ). 

Early-type galaxies in low-density environments ex- 
hibit small differences to those in clusters: at a given lu- 
minosity, the early-type galaxies in low-density regions are 
'-^1 — 3 Gyr younger and 0.1 — 0.2 dex more metal-rich 
than t hose in clusters ( Trager et al.l 2000al: iP oggianti et al.l 
| 200ll: IKuntschner et al.1 |2002|: ITerlevich fc Forbes 2002 



'Caldwell et al. 2003; Proctor et al 



, TT*^^ .2004al : IThomas et al.1 

2005.: Sanchez-BlazQuez et al.l |2006|) . In confiict with the 



above results, iGallazzi et al. I (|200d ). using a large sam- 
ple of early-type galaxies from the SDSS, found evidence 
that galaxies in low-density environments were less metal- 
rich than those in high-density environments. Intriguingly, 
there appears to be no e n vironmental depende nce of [a/Fe] 
IKuntschner et al.1 I2OO2I : IThomas et al] |2005| ). i ndicating 
that star formation in galaxies of a given mass occurs on 
the same timescale whether they are located in high-density 
or low-density environments. 

While it appears that the stellar populations of galaxies 
change from low-density environments to high-density envi- 
ronments, the density threshold at which this change occurs 
is uncertain. Studies of the star- formation rate of galaxies in 
and around clusters (|Lewis et al.l I2OO2I : IComez et all |2003| ) 
find an increase in the star-formation rate with increasing 
distance from the cluster centre, converging to the mean field 
rate at distances greater than ~ 3-Rvir. The critical projected 
density at which suppression of star-formation begins is un- 
certain but the environmental influences on galaxy proper- 
ties are believed not to be restricted to cluster cores, being 
effective in all groups where the density exceeds the critical 
value. The observed low rates of star formation well beyond 
the virialised cluster rule out physical processes associated 
with extreme environments (such as ram pressure stripping 
of disk gas) being completely responsible for the variations 
in galaxy properties with environment. 

This paper is the second in a series of papers aimed at 
studying the effects of environment and mass on the stel- 
lar populations of early-type galaxies. The first paper de- 
scribed our sample selection, observations, data reductions, 
and method of measuring line-strengths and estimating stel- 
lar population parameters within the framework of the Lick 
system. In this paper we present the results of this study 
utilising data from four clusters and their surrounds. The 
data extends to ~ 10i?vir, allowing us to probe the in- fall 
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regions of the clusters, which to date have been poorly stud- 
ied. The layout of the remainder of the paper is as follows. 
32]briefly describes the sample, observations and reductions, 
measurement of Lick indices and SPP estimation. Analysis 
of the absorption-line strengths is detailed in ij3]and the dis- 
tributions of the SPPs in §31 The correlations between the 
SPPs and velocity dispersion are discussed in ij5] while the 
contribution from these correlations to the intrinsic scatter 
in the parameter distributions are investigated in 3H1 We dis- 
cuss the results of this study in iJ7] and present a summary 
of the work in [JS] 

A Hubble parameter Ho ~ 70 km s~^ Mpc^^, matter 
density parameter Qm = 0.3 and dark energy density pa- 
rameter JIa = 0.7 are adopted throughout this work. 



2 FROM OBSERVATIONS TO STELLAR 
POPULATION PARAMETERS 

Details of the observations, data reductions, measurement of 
Lick indices and estimation of the corresponding SPPs are 
given in Harrison et al. (2010, submitted; hereafter referred 
to as Paper I). Therefore only brief descriptions will be given 
here. 

Observations of galaxies from four clusters (A930, 
A1139, A3558 and Coma) of varying richness and morphol- 
ogy were made with 2dF during the nights of 19-21 April 
2002. Observations of galaxies from the same four clusters 
and the structures around them were made with 6dF during 
the nights of 6-8 March 2003. The 6dF sample contained 
some galaxies in common with the 2dF sample, but mostly 
consisted of galaxies in the outer regions of each cluster. In 
two of the fields (A1139 and A930), the observations were 
offset from the cluster centre, allowing galaxies with cluster- 
centric distances of up to ~ 19 /i^q^ Mpc to be studied. The 
set-up of the two instruments can be found in Paper I. 

A brief discussion of the way galaxies were classified as 
early types or not is warranted here. This was done spectro- 
scopically with galaxies that showed signs of Ha (EW< —3.8 
A) or [OIII]A5007 A (EW< -0.4 A) emission being classified 
as star forming. We chose to classify a galaxy by its spectrum 
because comparisons to SSP models are (almost) meaning- 
less in galaxies with significant emission. In-fiU of the H/3 
absorption feature by nebular emission results in weaker 
H/3 line-strengths and incorrectly older ages . Some methods 
used to correct for this nebular in-fiU (e.g. lGonzale^ll993l : 
iTrager et al ]|2000a"^ have proved unsatisfactory, while oth- 
ers IjSarzi et a l. 2006) can only correct for relatively weak 
emission. For these reasons we eliminated all galaxies with 
signs of emission from our sample of early types. It must be 
noted that, in doing so, it is possible that we exclude from 
our stellar population analysis early-type galaxies that have 
had star formation triggered by the cluster environment. 

Other methods of classifying galaxies are not problem- 
free. Classifying by morphology is highly subjective and 
samples selected by colour-magnitude cuts still contain con- 
tamination by galaxies that would have been morphologi- 
cally classified as late types. More importantly, these meth- 
ods do not eliminate all emission-line galaxies, up to 30% 
of red sequence gala xies can show signs of LINER emission 
jCraves et al.ll2007l ). and so the interpretation of a stellar 
population analysis for these galaxies would be difficult. 



However, the absence of emission lines does not guar- 
antee that a galaxy has not had significant star formation 
in its recent past and so it is possible that our sample of 
early types contains a small number of post-starburst galax- 
ies. The intrinsic scatters in the line strength-a (Section 
13. 1|) and stellar population parameter-cr relations (Section 
I5.1|l are comparable to published results implying that if our 
sample contains recently star-forming galaxies then it is in 
no greater number than previous samples. 

The basic reduction steps, such as bias subtraction, 
spectrum extraction, flat-fielding, wavelength calibration, 
fibre throughput determination and correction, and sky 
subtraction, were performe d with the purpos e-built data 
reduction package s 2dfdr (|Colless et al.l I2001I ) and 6dfdr 



l|Jones et al. 



2004). Redsh ifts were measured using the pro- 



gram runz (ICoUess et al.ll2001^ while the IRAF task fxcor 



was used to measure velocity dispersions. 

Line-strengths were measured and transformed to the 



number of papers (e.g. Gonzaled 19931: 


Fisher et al. 19951: 


IWorthev & Ottavianilll997l: Traser et al. 


19981: Kuntschned 



(~ 9 A FWHM), and t he line-strengths we re measured us- 
ing the program indexf (jCardiel et al.|[l998l V 

A number of corrections were then applied to the mea- 
sured line-strengths to fully calibrate them to the Lick sys- 
tem. These corrections include a velocity dispersion correc- 
tion to account for the change in line-strength caused by 
velocity broadening, an aperture correction to account for 
the different linear sizes subtended by the different fibres 
at different redshifts, and, finally, applying any offsets that 
may arise due to the fact that the Lick spectra were not flux 
calibrated. This resulted in the measurement of the line- 
strength indices C24668, H/3, H/3g, [OIII]i, [0III]2, Fe5015, 
Mgi, Mg2, Mgb, Fe5270, Fe5335 and Fe5406 of the steUar 
populations within the inner ~ 1 kpc of each galaxy. 

The age, [Z/H], and [a/Fe] of a galaxy were obtained 
by comparing our mea sured line-strengths to the models 
of Irh omas et all (|2003l ) and accepting the c ombination of 
SPPs of the model with the smallest (jProctor et al.l 
l2004bl ). We interpolate the models logarithmically in steps 
of 0.02 dex in age and [Z/H], and 0.01 dex in [a/Fe]. Errors 
on the parameters were estimated b y using the consta nt 
boundaries as confldence limits (see lPress et al.lll992l ). 

In summary, we measured velocity dispersions, red- 
shifts, and line-strengths for a total of 416 galaxies: 158 
of these are early-type galaxies from Coma, A1139, A3558, 
or A930 that are located within the Abell radius (i.e. at a 
projected radial distance ^ 2 h^g Mpc) and comprise our 
cluster sample; 87 are early-type galaxies from the outskirts 
of these clusters at projected radial distances > 2 h^^ Mpc 
and comprise the cluster-outskirts sample; 168 galaxies were 
deemed to be star-forming, and make up the emission-line 
sample; the 3 remaining galaxies lacked the information nec- 
essary to classify them as either an early-type or emission- 
line galaxy. Of these 416 galaxies SPPs were estimated for 
219; 142 in the cluster sample and 77 in the cluster-outskirts 
sample. 
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Table 1. The line-strength-cr fits, Spearman rank correlation statistics, probabilities of a correlation, and the intrinsic scatter ((5intr; in 
mag) for the cluster galaxies. 



Index 


Zero-point 


Slope 


rs 


Probability 




C24668' 


-0.053 ±0.009 


0.067 ±0.004 


0.546 


1.000 


0.0096 ±0.0011 


H/3' 


0.167 ±0.009 


-0.044 ± 0.004 


-0.479 


1.000 


0.0075 ±0.0012 




0.238 ±0.011 


-0.054 ±0.005 


-0.442 


1.000 


0.0120 ±0.0022 


Fe5015' 


0.046 ± 0.008 


0.013 ±0.003 


0.284 


0.999 


0.0035 ± 0.0008 


Mgi 


-0.285 ±0.018 


0.180 ±0.008 


0.626 


1.000 


0.0198 ±0.0021 


Mg2 


-0.402 ±0.032 


0.299 ±0.014 


0.662 


1.000 


0.0200 ±0.0001 


Mgb' 


-0.130 ±0.012 


0.129 ±0.005 


0.675 


1.000 


0.0094 ± 0.0013 


Fe5270' 


0.044 ± 0.008 


0.016 ±0.004 


0.342 


0.999 


0.0045 ±0.0012 


Fe5335' 


0.008 ±0.011 


0.028 ± 0.005 


0.352 


1.000 


0.0081 ± 0.0009 


Fe5406' 


0.043 ±0.011 


0.010 ±0.005 


0.166 


0.921 


0.0035 ±0.0012 



3 ABSORPTION-LINE ANALYSIS 

3.1 Line-strength-cr relations in cluster early- type 
galaxies 

In early-type galaxies the strengths of certain absorption 
features are observed to be correlated with velocity disper- 
sion. Such correlations are important as they provide links 
between a galaxy's dynamical and chemical evolution; i.e. 
between galaxy mass, metallicity, and abundance ratios. 

Historically, the Mg-a relation was studied using Mg2, 
which is measur ed in magnitudes of absorbed flux. More re - 
cent studies (e.g. lColless et aO 19991 : iKuntschner et al ] |200ll ) 
have used Mg6, which is a more rehable index, as it is nar- 
rower, but is measured as an equivalent width in A. These 
studies converted indices measure d in equivalent wid ths to 
magnitudes (for the conversion see IColless et ai][l999l l. pri- 
marily to allow comparison with older studies that used 
Mg2. We continue this practice and indices that have been 
converted to magnitudes will be denoted with a prime sym- 
bol (e.g. Mg6')- 

The variations of the index line-strengths with veloc- 
ity dispersion for the cluster sample are shown in Figure [1] 
These plots contain the combined data from all four clus- 
ters, as Kolmogorov-Smirnov (KS) 2D two-sample tests in- 
dicate consistency of the line-strength-cr distributions be- 
tween clusters, in that most distributions differed only at 
the l-2cr level. That these relations are found to be con- 
sistent in the four clusters, which span the ranges of Abell 
richness classes and Bautz-Morgan (EM) morphologies, is 
remarkable. This suggests that the cluster environment has 
little effect o n the l ine-strengths of early-type galaxies (see 
ICoUess et all (|l999l ) for a study of the Mg-cr relation as a 
function of cluster environment). Table [1] shows the results 
of linear fits to the data, accounting for errors in both quan- 
tities, along with the Spearman rank correlation statistic 
(rs), the probability that the two quantities are correlated 
and the estimated intrinsic scatter about each relation. 

We assume the intrinsic scatter about these relations 
is Gaussian and estimate it by maximising the logarithmic 
likelihood 



log ^ = X! 



: exp 



2S? 



(1) 



where the Ali axe the residuals to the fit for index / and 
Si = Sif + {ai 5 log CTi)^ + 5,?„tr) where Sit and S log ai are the 



individual errors on the index / and log a, ai is the slope of 
the relation for index /, and (5intr is the estimated intrinsic 
scatter. These estimated intrinsic scatters are listed in Table 

m 

With the exception of Fe5406, correlations are found 
at > 3cr confidence level for all indices. Strong correlations 
(rs > 0.5) are found for C24668, Mgi, Mg2 and Mgb, while 
a weak correlation (rs < 0.3) is found for Fe5015. Except 
for H/3 and H/3g, which are moderately anti-correlated, all 
other indices are moderately correlated. The slope of the fit 
to the Fe5406 data is very flat and consistent with being 
zero. H/3 and H/3g are the only indices used here that are 
more sensitive to age effects and the only ones that exhibit 
an anti-correlation; all the other indices are more sensitive 
to [Z/H] effects and either exhibit a positive correlation or 
no correlation. 

Young stellar populations exhibit strong Balmer ab- 
sorption features, which weaken with age. Therefore, a sim- 
ple interpretation of the anti-correlation of H/3 and H/3g 
with velocity dispersion is that more massive galaxies are 
older. The simple interpretation of the positive correlation 
between the metal-sensitive indices and velocity dispersion 
is that more massive galaxies are more metal-rich. The usual 
explanation of this trend is that massive galaxies are better 
able to retain their heavy elements due to their larger gravi- 
tational potential, a scenario that arises na turally in galactic 
wind models fe.g. lArimoto fc Yoshii|[l987h . Alternatively, a 
vari able IMF could l ead to a similar mass-[Z/H] relation. 
iKoeppen et al.l (|2007l ) have developed a model where the ef- 
fective upper mass limit of stars is lower in galaxies with a 
low star-formation rate. This reduces the number of SNII, 
and hence the [Z/H], in low mass galaxies and leads to a 
similar mass- [Z/H] relation. 

We compare our line-strength-cr relations to vari- 
ous estimates from the literature in Table [2] These esti- 
mates are based on a sample of early-type galaxies from 
the Fornax cluster (jKuntschneil 2000), a sa mple of red- 
seque nce galaxies from numerous low-z clusters l|Nelan et al.l 
|2005| ). a magnitude- limited sample of early-ty pe galax- 
ies drawn from the SDSS (jClemens et al.l l2006h . a sam- 
ple of oarly-typo galaxies in high-density environments 
( Sanchez-Blazquez ct al. 2006), a magnitude-limited sam- 
le of early-type g alaxies in high-density environments 
Ogando et al.l [20081 ) . and a samp le of early-type galax ies 
from the core of the Coma cluster (|Matkovic et al.ll2009l '). 
Generally, our results compare well with those in the 
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Figure 1. The variations (in magnitudes) of the index hne-strengths with velocity dispersion for the cluster sample. Most indices exhibit 
a tight relation with strong correlations found for C24668, Mgi, Mg2 and Mg6, while H/3 and H/3q exhibit moderate anti-correlations. 



literature. We find that in the majority of cases tfie slopes 
agree at tfie 2a level or better. Of those comparisons that 
disagree by more than 3a, we note that almo s t half are com- 
parisons to the estimates of IClemens et all (|2006l ) and the 
rest consist almost entirely of comparisons of Mgi and Mg2 . 

The discrepancy with I Clemens et all (j2006l ) possibly 
arises due to differences in the aperture corrections applied 
to each dataset. We correct our line-strength measurements 
to fixed physical size of ~ 1 kpc, independent of the ve- 
locity dispersio n of the galaxy. The correction applied by 
IClemens et al.l . however, is to a physical size of ro/10 and 



is a function of velocity dispersion. Confusing the issue is 
the fact that some of their slopes are steeper than ours (e.g. 
Fe5406), some are shallower (e.g. C24668) and some are in 
the opposite sense (e.g. Fe5015). 

The discrepancy between our estimated slopes of Mgi 
and Mg2 and those in the literature is a little easier to under- 
stand. Both these indices have very broad definitions and so 
their line-strengths are highly-sensitive to the shape of the 
continuum. The fact that we divided out our continuum is 
probably the reason for this discrepancy. This is one of the 
reasons why we excluded these two indices from our fits to 
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Table 2. A comparison of the line-strength— cr relations to those found in the literature. 



Reference 




C24668' 


H/3' 




Fe5015 








Zero-point 


Slope 


Zero- point 


oiope 


Zero-point 


oiope 


This study 






n nfi7+n 004 

\J.\J\J 1 _I_U . VJU'rl 


0.167±0.009 - 


-0.044±0.004 


0.046±0.008 


0.013±0.003 


Kuntschner (2000) 




_n 1 1 0-1-0 049 


nQD-i-n 01 s 


0.106±0.015 - 


-0.020±0.007 


0.002±0.019 


0.036±0.008 


Nelan et al. ('2005) 






f)7'i-l-n 002 




-0.041±0.001 




0.015±0.001 


Clemens et al. (2006) 




0.044±0.003 


0.025±0.003 


0.227±0.002 - 


-0.071±0.002 


0.127±0.002 - 


0.023±0.002 


Sanchez-BlazQuez et al. 


(2006) 


-0.071±0.018 


0.067±0.008 


0.087±0.009 - 


-0.012±0.004 


0.025±0.011 


0.021±0.005 


Ogando et al. (2008) 








0.156±0.013 - 


-0.042±0.006 


0.053±0.010 


0.011±0.004 


iviaLKOvic cL ai. i zuuy j 




-0.044±0.014 


0.060±0.014 


0.116±0.009 - 


-0.024±0.009 


0.060±0.007 


0.005±0.007 


Reference 




Mgi 




Mg2 




Mg6' 








Zero-point 


Slope 


Zero-point 


Slope 


Zero-point 


Slope 


This study 




-0.285±0.018 


0.180±0.008 


-0.402±0.032 


0.299±0.014 


-0.130±0.012 


0.129±0.005 


Kuntschner (2000) 




-0.158±0.035 


0.136±0.015 


-0.127±0.054 


0.191±0.023 


-0.056±0.044 


0.102±0.020 


Nelan et al. (2005) 






0.121±0.003 




0.189±0.003 




0.134±0.002 


Clemens ct al. (2006) 




-0.245±0.010 


0.163±0.004 


-0.198±0.013 


0.208±0.006 


-0.155±0.007 


0.142±0.007 


Sanchcz-BlazQucz et al. 


(2006) 










-0.050±0.019 


0.091±0.008 


Oeando et al. (2008) 




-0.159±0.020 


0.131±0.009 


-0.153±0.029 


0.194±0.013 


-0.095±0.018 


0.112±0.008 


Matkovic et al. (2009) 




-0.204±0.036 


0.141±0.017 


-0.143±0.045 


0.178±0.020 


-0.023±0.015 


0.080±0.015 


Reference 




Fe5270' 


Fe5335' 


Fe5406' 






Zero-point 


Slope 


Zero-point 


Slope 


Zero-point 


Slope 


This study 




0.044±0.008 


0.016±0.004 


0.008±0.011 


0.028±0.005 


0.043±0.011 


0.010±0.005 


Kuntschner (2000) 




0.024±0.020 


0.029±0.009 


-0.017±0.020 


0.043±0.009 


0.023±0.026 


0.023±0.012 


Nelan et al. (2005) 






0.017±0.001 




0.023±0.001 




0.018±0.001 


Clemens et al. (2006) 




0.085±0.003 - 


0.001±0.003 


0.014±0.003 


0.027±0.003 


-0.002±0.003 


0.031±0.003 


Sanchez-BlazQuez et al. 


(2006) 


0.027±0.007 


0.024±0.003 


-0.001±0.013 


0.035±0.006 






Osando et al. (2008) 




0.049±0.008 


0.015±0.004 


0.029±0.009 


0.020±0.004 


0.034±0.009 


0.018±0.004 


Matkovic et al. (2009) 




0.092±0.009 - 


0.007±0.009 


0.069±0.008 


0.001±0.008 







the models in estimating the SPPs (see Section 6 in Paper 
I). 

Most comparisons with lMatkovic et al.l ()2009l ) agree at 
better than the 3(7 level, the exceptions being Mg2 and Mgb, 
but since both samples include early-types from the core of 
the Coma cluster it is worthwhile to make a direct compar- 
ison of these galaxies alone. Overall, this results in better 
agreement between the slopes, although the change is al- 
most negligible in most cases. Three indices (C24668, Fe5270 
and Fe5335) show worse agreement but the largest change 
in significance is only 0.24cr. All other indices show better 
agreement; again, most changes are small except in the case 
of H/3 where the significance changes by 1.38cr from 2.03(t 
to 0.65cr. 

The intrinsic scatters found for our line-strength- 
a relations are given in Table [1] Most o f our esti- 
mates agree well with those of iKuntschneil (|2000t ) and 



ISanch cz-Blazguc z et al.l (|2006l ). with perhaps the exceptions 
of H/3, Fe5015, Fe5270 and Fe5335. We find an intrinsic scat - 
ter in H/3 of 0.008 mag ( as does Sanchez-Blazguez et al.h . 
which is double that of IKuntschneil . For the three Fe 
indices, our scatters are larger than those found by 
ISanchez-Blazguez et al.l who find negligible-to-no scatter in 
these relations. Our esti mates for th ese indices are in good 
agreement with that of IKuntschneil except for Fe5335 for 
which we find almost double the scatter of that found by 
this author (0.008 mag compared to 0.005 mag). 

Given that the galaxies from this study are drawn 
from four different clusters, with varying richness classes 
and BM classifications, the small intrinsic scatters about 



the line-strength-o" relations is truly remarkable, and im- 
plies that the stellar populations in cluster early-type 
galaxies are homog eneous. This result agrees with that of 
IColless et al.l l|l999l ). who found no correlation between the 
Mg-o" zero-point and cluster velocity d ispersion. X-ray lumi- 
nosity, or X-ray temperature (see also Bernard! et al.l 19981 : 
IWorthev fc CollobertI 12003': 'Sanc hez-Blazguez et al.ll2006l ). 
For gala:xies in clusters, Jorgcnsc n et al. I (|l996l ) found that 
the Mg-(T residuals are weakly correlated with local density. 
However, since the residuals are correlated with cluster ve- 
locity dispersion and not projected cluster-centric distance, 
the interpretation is of a correlation between the Mg-cr zero- 
point and cluster mass, albeit a weak one. 



3.2 Line-strengths— «r relations in the 
cluster-outskirts 

Using the wide field of 6dF and, in the cases of A1139 and 
A930, offsetting the field centres with respect to the cluster 
centres allowed us to obtain data on galaxies well outside 
the clusters. We obtained line-strengths of galaxies out to 
a projected cluster-centric radius of ~ 19 K^q Mpc, which 
corresponds to ~ 107?vir. We note that the average distance 
between Abell clusters is ~ 30 hj^ Mpc. 

The trends with velocity dispersion of the line-strengths 
in the cluster-outskirts sample (i.e. those galaxies outside 
^?Abeii) are shown in Figure[2l The solid line in each panel is 
our linear fit to the data (accounting for errors in both quan- 
tities) while the dashed line is the fit to the cluster galaxies. 
In Table |3] we list the details of the fits, the Spearman rank 
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Figure 2. The variations (in magnitudes) of the index hne-strengths with velocity dispersion for the cluster-outskirts sample. The solid 
line is the linear fit to the data while the dashed red line is the fit to the cluster sample. The slopes of all relations are consistent with 
those found in the cluster sample. 



correlation statistics, the probabilities of a correlation and 
the estimated intrinsic scatters about the relations. The last 
column in this table shows the number of standard devia- 
tions by which the slopes of the cluster-outskirts sample re- 
lations differ from the cluster sample relations; all indices are 
found to have slopes consistent with those found for galax- 
ies in the cluster sample. However, the only indices that 
show correlations significant at the 3cr level in the cluster- 
outskirts sample are Mg6 and H/3g, and while the strength 
of the correlation remains the same for H/3g it is reduced 



for Mg6 from a strong correlation in the cluster sample to a 
moderate correlation here. 

The change in the Fe indices is interesting. Fe5015 
changes from being weakly correlated in the cluster sample 
to uncorrelated in the cluster-outskirts sample. The other 
three Fe indices remain moderately-to- weakly correlated (al- 
though with reduced significance) but their slopes are all 
consistent with being zero. In the cluster sample only Fe5406 
had a slope that was consistent with being zero. In addition 
all four Fe indices have decreased intrinsic scatters about 
their fits, as does Mgi. C24668, II/3, H/3g and Mg6 all show 
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Table 3. The line-strength— cr fits, Spearman rank correlation statistics (rg), probabilities of a correlation, and intrinsic scatters {5intr; in 
mag) for galaxies in the cluster-outskirts sample. The number of standard deviations by which the slopes of the cluster-outskirts sample 
relations differ from the cluster sample relations ((Tdiff ) is given in the last column. 



Index 


Zero-point 


Slope 


rs 


Probability 


'Jintr 




C24668' 


-0.064 ±0.017 


0.069 ± 0.007 


0.260 


0.985 


0.0200 ± 0.0020 


0.3 


H/3' 


0.171 ±0.017 


-0.044 ± 0.007 


-0.275 


0.988 


0.0184 ± 0.0045 


0.0 




0.264 ±0.021 


-0.066 ±0.009 


-0.322 


0.997 


0.0200 ± 0.0003 


1.2 


Fe5015' 


0.053 ±0.016 


0.010 ±0.007 


0.086 


0.413 


0.0029 ±0.0018 


0.4 


Mgi 


-0.316 ±0.056 


0.200 ±0.024 


0.748 


0.995 


0.0142 ±0.0038 


0.8 


Mg2 


-0.363 ±0.078 


0.278 ±0.034 


0.650 


0.978 


0.0200 ± 0.0036 


0.6 


Mg6' 


-0.108 ±0.023 


0.115 ±0.010 


0.392 


1.000 


0.0200 ±0.0012 


1.3 


Fe5270' 


0.016 ±0.035 


0.027 ±0.015 


0.410 


0.814 


0.0000 ±0.0011 


0.7 


Fe5335' 


0.003 ± 0.030 


0.028 ±0.013 


0.548 


0.990 


0.0000 ± 0.0020 


0.0 


Fe5406' 


0.024 ± 0.030 


0.017 ±0.013 


0.405 


0.939 


0.0027 ±0.0017 


0.5 



an increase in the intrinsic scatter, while for Mg2 it remains 
the same. 

There appears to be a large number of galaxies that 
have low index strengths compared to those expected from 
the relations found in the cluster sample. For Fe5015 the 
spread is consistent with what is found in the cluster sam- 
ple. The galaxies with low H/3 and H/3g strengths, which 
also are found in the cluster sample but in smaller numbers, 
are possibly those that suffer from nebular in-fiU. Most of 
the galaxies in the cluster-outskirts sample were observed 
with 6dF and so have no details of Ha emission because 
the wavelength coverage was insufficient. Therefore, it is not 
unexpected that the level of this possible contamination by 
star-forming galaxies in this sample is slightly higher than 
in the cluster sample. 

In the cases of C2 4668' and Mgb, there does appear to 
be a real sub-population of galaxies that are genuinely offset 
from the relations found in the cluster sample, having weaker 
line-strengths for a given velocity dispersion. This is possible 
if galaxies in clusters find it easier to retain their heavier 
elements than those in their outskirts, due to the dense intra- 
cluster medium retarding the development of galactic winds 
and the expulsion of the heavier elements. 

Another possibility is that the weakened line-strengths 
are caused by the aperture correction. Early type galax- 
ies exhibit strong metallicity and age gradients (e.g. 
iDavies et al.|[l993i ) and indices measured in massive galax- 
ies will sample a smaller fraction of the effective radius than 
those in less massive galaxies. If lower mass galaxies are 
found preferentially in the cluster outskirts then this could 
lead to a change in the relations in the two regions. These 
aperture corrections are necessary to allow a fair compar- 
ison between galaxies in clusters at different redshifts and 
between galaxies in the same cluster but observed with fibres 
that subtend different angles (2dF~ 2.1" and 6dF~ 6.7"). 
However, we do not find that lower mass galaxies are prefer- 
entially found in the cluster-outskirts sample. In fact, while 
the distributions of magnitudes are similar in shape in both 
regions (because we deliberately targeted brighter galaxies), 
the peak of the distribution in the core is shifted to fainter 
magnitudes relative to the outskirts sample by ~ 1 mag, 
i.e. on average the cluster-outskirts galaxies are brighter. 
Therefore this can not be the reason for the weakening of 
the relations. 



3.3 The radial distribution of star-forming 
galaxies in clusters 

In the analyses so far, the included galaxies were limited to 
those that showed no significant sign of on-going star for- 
mation, as determined from Ha or [OIIIJA5007A emissions. 
This allows us to be reasonably confident that their H/3 
absorption feature is free from in-filling caused by nebular 
emissions. However, from these previously excluded galax- 
ies we can determine where in the cluster and its surrounds 
star-forming galaxies reside. 

In the left panel of Figure [3] we show, for our entire 
sample of galaxies, the H/3 line-strengths as a function of 
projected cluster-centric distance. Circles represent galax- 
ies that are classified as early-types, while crosses repre- 
sent emission-line galaxies (and were up to this point ex- 
cluded from our analysis). The distribution of H/3 line- 
strengths forms a ridge at ~ 1.8 A. Most of the galax- 
ies with H/3 line-strengths weaker than ~ 1.4 A are those 
classified as emission-line gala^des, although a few galax- 
ies classified as early-types also lie below this ridge line. 
It is possible that these galaxies have weak H/3 due to in- 
filling and yet show no significant sign of [OIII]A5007A emis- 
sion. This confirms that while correcting for H/3 e mission by 
the strength of the [OII I]A5007A emission line jGonzalej 
ll993l : lTrage~ et al. I l2000al lbl) may be acceptable in a statisti- 
cal sense, it is not reliably applicable to individual galaxies, 
which may show moderate to strong H/j emissio n and little 
or no [OIII]A5007A emission iNelanetal .1120051 ). 

Figure [3] shows an increasing scatter in H/3 line- 
strengths with increasing radius. Not only is there a scatter 
to lower values (an indication of older ages or possibly H/3 
in-filling in galaxies that have experienced recent star for- 
mation), there is also a scatter to higher values, and thus 
younger ages, for early-type galaxies. The red line in the left 
panel of Figure |3] shows how the mean H/3 line-strength (cal- 
culated for galaxies within a 0.4 dex radius bin at steps of 
0.05 dex) changes with radius. We see that within a radius 
of 2 hjQ Mpc (i.e. ~ /?Abcii) the mean value is a constant 
(~ 1.8 A), while outside this radius it steadily decreases to a 
value of ^ 1.4 A. Sin ce the oldest isochrone (15 Gyr) in the 
iThomas et"all (|2003l ) models has H/3 strengths ~ 1.6-2 A we 
conclude that these low line-strengths, for the galaxies clas- 
sified as early types, represent nebular in-fiU in recently star- 
forming galaxies, which are preferentially found in the outer 
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Figure 3. Left: The variation in H/3 line-strength as a function of projected cluster-centric radius. Circles represent early-type galaxies 
and crosses represent galaxies that are classified as emission-line galaxies. The red line shows how the mean H/3 line-strength changes 
with radius and indicates that on average galaxies in the outskirts of clusters are ^ 0.4 A weaker in H/3 than those in the centre of 
clusters. The blue line shows how the H/3 line-strength evolves in a galaxy as it ages after an instantaneous burst of star formation (the 
direction is given by the blue arrow; see text for details). Right: The fraction of galaxies with H/3 1.4 A as a function of projected 
cluster-centric radius. 



regions of the cluster environment. This conclusion is consis- 
tent with a prev ious study of galax y star-formation rate and 
environment bv lLewis et al.l ( 20021 ) . which found increasing 
star-formation rates with increasing distance from the clus- 
ter centre that conver ge to field rates at distance s greater 
than ~ Siivir (see also iGomez et al.l[200^ : iXauffmann et al.l 
|2004 ). 

The distribution of galaxies in this figure can be thought 
of as tracking the evolution of a star-forming galaxy. A 
galaxy falling into the cluster undergoes a burst of star for- 
mation, which shifts the galaxy to H/3 emission (i.e. negative 
line-strengths). As the episode of star formation progresses 
and the galaxy moves deeper into the cluster, the galaxy 
gradually shifts to stronger H/3 absorption (positive line- 
strengths) . Then, once star formation has ceased, the galaxy 
settles back onto the ridge. This evolution is represented 
schematically by the blue line in the figure. The dotted seg- 
ment of the line represents the galaxy as it moves from be- 
ing a star-forming galaxy to a post-starburst galaxy (the H/3 
emission line-strengths are indicative only). The solid seg- 
ment of the line shows how, according to the iThomas et al.l 
models, the H/3 line-strength changes in a galaxy with 
[Z/H]=0.35 dex and [a/Fe]=0.2 dex as it ages from 1 Gyr to 
15 Gyr. The mapping between projected cluster-centric ra- 
di us and time is arbit r ary. T his figure is similar to Figure 10 
in I Couch fc Sharpie^ (|l987l ) , which shows the evolutionary 
track of a star-bursting galaxy in H5-colour space. 

The right panel of Figure [3] shows the fraction of galax- 
ies (in the same bins as used for calculating the mean in the 
left panel) that have line-strengths less than 1.4 A. In the 
core of the cluster none of the galaxies fall below this value. 



while at large radii the fraction increases rapidly to the limit 
of our data, where ~ 40% of galaxies have line-strengths 
less than this value. We conclude that the cluster core is 
relatively free from young galaxies and galaxies that have 
experienced recent star formation, and that these galaxies 
are found more commonly outside i?Abeii. 



4 STELLAR POPULATION PARAMETER 
DISTRIBUTIONS 

4.1 The cluster distributions 

The distributions of SPPs from the four clusters combined 
are shown in Figured Marginal distributions are plotted for 
each of the parameters and the median values are marked 
as dotted lines. Galaxies from Coma are shown as blue cir- 
cles, from A1139 as green squares, from A3558 as yellow 
triangles and from A930 as red crosses. Errors are shown 
for individual galaxies but two points must be kept in mind. 
Firstly, the age and [Z/H] errors are correlated and so the 
error bars here do not accurately reflect the true shape of 
the confldence contours. For the two other combinations of 
parameters the errors are much less correlated and so the er- 
ror bars are a good representation of the confidence contours 
(see Paper I). Secondly, a galaxy's error estimate cannot be 
larger than its distance to the edge of the model grid. This 
is only an issue for age estimates since we quote two-sided 
errors. So, if a galaxy's error bar reaches 15 Gyr it should 
be considered a lower limit only. 

Looking at the distributions as a whole, we note that 
very few galaxies have [Z/H] less than solar with most having 
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[Z/H] Age (Gyr) [a/Fe] 

Figure 4. The distributions of the SPPs for all galaxies in the cluster sample. Marginal distributions are shown for each parameter and 
the dotted lines show the median values of each parameter. Galaxies from Coma are shown as blue circles, from A1139 as green squares, 
from A3558 as yellow triangles and from A930 as red crosses. 



Table 4. The median SPP values for each cluster and the number 
of galaxies for which the parameters were measured. 



Cluster 


[Z/H] 


Age 


[a/Fe] 




Coma 


0.25 ±0.02 


6.6 ±0.6 


0.27 ±0.01 


65 


A1139 


0.19 ±0.08 


4.9 ±1.6 


0.32 ± 0.04 


8 


A3558 


0.29 ±0.04 


5.7 ±0.9 


0.29 ± 0.02 


61 


A930 


0.20 ±0.08 


6.4 ±2.1 


0.27 ±0.05 


8 



0.1 < [Z/H] < 0.4 dex. Similarly, [a/Fe] is mostly greater 
than solar and lies in the range 0.2 < [a/Fe] < 0.4 dex. The 
bulk of the galaxies have old ages and almost all lie in the 
range 4 < age < 15 Gyr. The distributions we find for the 
cluster sample are in general agreement with those found 
by previous authors (e.g. [Gonzalez 1993: Jorgensen 1999. 
Kuntschneill2000l:lTrager et al.ll2000briPoggianti et al.ll200 



Thomas et al.ll2005l : ICoUobert et al.ll200i 

Looking at the SPP distributions of each cluster in- 
dividually we find few differences. The median values and 
errors for the SPPs in each cluster are given in Table ID 
The errors on the median values were estimated from Monte 
Carlo bootstrap simulations in the following manner. For 
each galaxy in a cluster, we randomly draw index values 
from the galaxy's index error distributions, using the same 
set of indices that were used to estimate the parameters 
originally. These index values are then converted to param- 
eter estimates and the median values are determined. This 
process is repeated 10,000 times and the RMS errors on the 
median values are computed. 

The distributions of SPPs from cluster-to-cluster show 
remarkable similarity, with no significant differences between 
the median values of any of the four clusters. However, to 
conclusively detect any differences in the age and [Z/H] dis- 
tributions, a two-sample 2D KS test is required, due to the 
correlated errors between age and [Z/H]. For the [a/Fe] dis- 
tributions a two-sample ID KS test is sufficient. The results 
of these tests are that the joint distributions of age and 



[Z/H], and the distributions of [a/Fe] in each of the four 
clusters are consistent with each other. 

Looking at the combined distributions of the cluster 
sample (142 galaxies in total), we find that the median [Z/H] 
is 0.25 ± 0.14 dex, the median age is 6.01^2 g Gyr, and the 
median [a/Fe] is 0.28±0.11 dex. The errors on these median 
values were determined the same way as for the clusters 
individually. These median values are shown as dotted lines 
in the marginal distributions in Figure [l] 

There are some galaxies that have line-strengths that 
fall outside the ranges predicted by the stellar population 
models, in that there are six galaxies with age=15 Gyr, four 
with [Z/H]=0.59 dex, and two with [a/Fe]=-0.19 dex. Two 
galaxies have more than one parameter on the edge of the 
models and so in total there are 10 such galaxies. These 
galaxies are assigned the combination of SPPs that most 
nearly fits the combination of line-strengths, as mentioned 
in Paper I. It is hard to determine the origin of this differ- 
ence between the data and the models, which could be due 
to errors in the stellar population models (it should be noted 
that the models provide no estimate of systematic uncertain- 
ties), or due to errors in the observations or data reduction; 
only a small percentage of galaxies are affected, with ~ 90% 
of galaxies having self-consistent model parameter fits. 

For the galaxies with age=15 Gyr, it is possible that 
some of these galaxies are affected by nebular H/3 emis- 
sion, despite the fact that care was taken in eliminating such 
galaxies. As was noted in Paper I, it is possible for galax- 
ies to have FI/3 in emission and no detectable [OIII]A5007A 
(|Nelan et al.ll2005l ). If this H/3 emission is undetectable, be- 
cause it is swamped by H/3 absorption, and there is no 
[OIII]A5007A emission, then such a galaxy will not be clas- 
sified as an emission-line galaxy and will not be excluded 
from the cluster sample. 

It is also evident that there appears to be an anti- 
correlation between [Z/H] and age, in the sense that younger 
galaxies are more metal-rich. This might in principle be due 
to the non-orthogonal nature of the [Z/H]-age grids pro- 
duced by the stellar population models, which mean that 
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Figure 5. The distributions of the SPPs, estimated from H/3, Mgfe and Fe5335 only, for all galaxies in the cluster sample. The dotted 
lines show the median values of each parameter. 



the errors on these two quantities are correlated. The effect 
is that an increase in a [Z/H] -sensitive index results not only 
in a increase in [Z/H] but also in a decrease in age; corre- 
spondingly, an increase in an age-sensitive index results in a 
decrease in age and an increase in [Z/H]. Due to the irreg- 
ular shape of the grid, this effect varies with location over 
the grid. 

That a degree of anti-correlation between age and [Z/H] 
is introduced by the correlated errors is not disputed; how- 
ever, whether this ac counts for the anti- c orrelation entirely 
is still contentious (IColless et al.l Il999l: J0rgensenl Il999al: 
Trager et al. 2000al: Poggianti et al.l 2001 : Kuntschner et al.l 



200li: iTerlevich fc ForbesI l2002l: iProctor fc Sansoml 1200^ 



Bernardi et all l2005l : ISanchez-Blazguez et al.l 120061 '). Some 
studies find that there exists a moderate anti-correlation 
over and a b ove that introduced by the errors (e.g. 
IColless et all Il999l ) while others find that no correla- 
tion exists once the e rrors are taken into account (e.g. 
iKuntschner et al] 1200 lI V To determine whether the anti- 
correlation between age and [Z/H] is real or not, we perform 
a simple test which involves comparing the degree of anti- 
correlation in the observed joint age-[Z/H] distribution with 
those obtained by Monte Carlo simulations. 

Due to the need for an automated process to convert 
SPPs to line-strengths (and vice versa) it is not feasible 
to use the method of SPP estimation that makes use of 
all available indices. Therefore, for this exercise, we rely on 
only three indices H/3, Mg6 and Fe5335. Combinations of 
these three indices are commonly used in line-diagnostic di- 
agrams to estimate SSPs. With this combination of indices, 
the age distribution is approximately exponential and both 
the [Z/H] and [a/Fe] distributions are approximately Gaus- 
sian (see Figure [S}. 

We start by generating a sample of galaxies with age, 
[Z/H], and [a/Fe] estimates (parameter triples); the [Z/H] 
and [a/Fe] were drawn from the distributions defined by 
the medians of the observed values and their intrinsic scat- 
ters and the age was drawn from an exponential distribution 
(see Section |6] for details on how the intrinsic scatters and 
the e-folding of the exponential age distribution were de- 
termined). These parameter triples are then converted to 
H/3, Mg6, and Fe5335 index values (index triples), which 
are then perturbed using the observed indices error distri- 



butions. These perturbed index triples are then converted 
back to parameter triples and the Spearman rank correlation 
statistic rs is calculated for age and [Z/H]. This procedure 
is carried out 10,000 times. From these Monte Carlo simu- 
lations we determine the probability of obtaining by chance 
a correlation greater than that found for the observed data, 
and we find that there exists a real anti-correlation between 
age and [Z/H] (over and above the error-induced correlation) 
that is significant at the > 3cr level. 



4.2 The distributions in the cluster-outskirts 

The distributions of SPPs in the cluster-outskirts sample 
(red crosses) are compared to those from the clusters (black 
circles) in Figure [5] For the sake of clarity, only the errors 
for the cluster-outskirts sample are shown. The dashed lines 
in the marginal distributions are the median galaxy parame- 
ters for the cluster-outskirts sample and the dotted lines are 
for the cluster sample. The gala:xies in the cluster-outskirts 
sample have a median [Z/H] of 0.21 ±0.27 dex, a median age 
of 6.61^:6 Gyr, and median [a/Fe] of 0.27 ± 0.21 dex. Com- 
pared to the combined cluster sample we detect no difference 
in the median values for the galaxy parameters, although the 
errors are substantial. 

Looking at the distributions we see that there is very 
little difference between the galaxies in this sample and those 
in the clusters, an assessment that is confirmed by KS testing 
(a two-sample 2D KS test in the case of the age and [Z/H] 
distributions and a two-sample ID KS test for the [a/Fe] 
distributions) . 

The fact that the [a/Fe] distributions are the same in 
the clusters and their outskirts is intriguing, given that the 
cores of galaxy clusters are thought to result from regions of 
high ove r-density where the process o f star formation occurs 
rapidly jK auffman n fc Charlotlll998l: ISchindler et al.ll2005l : 
iRomeo et al.,2005.;iDe Lucia et al.ll2006l ). Our result is. how- 
ever, in agreement with previous s tudies on the [a/Fe] in 
low- and h igh-density environments (jKuntschner et cil]|2002l : 
lThomas~et al. 2005). These results suggest that, in all envi- 
ronments, elliptical galaxies form on similar timescales. 

The often-found differences in low- and high- 



densit y environments b etwee n the age llTrag:er et af] 
l2000al : iPoggianti et~al] I2OOII : IKuntschner et al.l l200a 
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Figure 6. Comparison of the SPPs in the cluster-outskirts sample (black circles) to those in the cluster sample (red crosses). The 
marginal distributions for the cluster sample are shown as open histograms and for the cluster-outskirts sample as black histograms. The 
dashed line in the marginal distributions represents the median parameter values in the cluster-outskirts sample and the dotted line that 
of the cluster sample. 



Tcrlcvic h fc Forbe ^'2002': Caldw ell et al ]l2003';P roct or et al 
2004a ; Thomas et al.. .2005 ) and [Z/H] (|Kuntscliner et al 



20021 : iThomas et allbOOsl T with galaxies on average being 



^ 2 Gyr younger and more metal-rich in lower-density 
environments, are not reproduced here. Although, given the 
size of our errors it is not surprising that we do not detect 
such a small difference. These results are usually found 
comparing truly isolated galaxies with cluster galaxies. As 
the galaxies in the cluster-outskirts sample are actually 
drawn from the outer regions of the clusters and the 
structures surrounding them (out to a projected radial 
distance of ~ 19 h^^ Mpc or ~ lOJJvir), it is possible 
that the contrast between the average densities of each 
environment is not sufficient to reveal any differences in the 
stellar populations, and that the change occurs at a lower 
density threshold (e.g. between isolated galaxies and those 
in groups/filaments). 

Figure [7] shows the variations of the SPPs with pro- 
jected cluster-centric distance, normalised to the cluster 
virial radius (i?vir). The lines in these plots show the mean 
values in bins containing 15 galaxies extending out to al- 
most 10-Rvir. The error bars show the standard error of the 
mean within each bin. We find no evidence of any signif- 
icant trends, although it does appear as if [a/Fe] and age 
increase moving towards the centre of the cluster, imply- 
ing that the galaxies in the cluster cores are older and have 
shorter star-formation timescales than those in the cluster 
outskirts. Despite the lack of clear trends there are two other 
interesting aspects of these plots. Firstly, there appears to be 
a decrease in [Z/H] outside ~ 27?vir. Secondly, at a distance 
slightly less than the cluster virial radius there appears to be 
a dip in the average age. Such a dip would be expected in a 
scenario where interaction with the ICM triggers star forma- 
tion in an in-falling galaxy, resulting in a stellar population 
characterised by a young age and extended star-formation 
timescale. This decrease in age begins at a radius > 3J?vir 
suggesting that the inffuence of the cluster extends to large 
distances. We also note that all galaxies with ages < 2 Gyr 
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Figure 7. The trend of mean age (top), mean [Z/H] (middle), 
and mean [a/Fe] (bottom) as a function of cluster virial radius. 
The mean is calculated in bins containing 15 galaxies. 



are located near i?vir and the age of the youngest stellar 
population at a given radius increases moving towards the 
cluster core. 



5 THE PARAMETER-o- RELATIONS 
5.1 The cluster relations 

The variations of the SPPs with velocity dispersion are 
shown in Figure [8] All parameters are found to be mod- 
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log(CT [km s ']) log(t7 [km s ']) log(CT [km s ']) 



Figure 8. The variations of the SPPs log(age) (left panel), [Z/H] (middle panel) and [«/Fe] (right panel) with log it of the galaxies in 
the cluster sample. All three parameters are moderately correlated with log a. The red line in each panel is the linear fit to the data. 



erately correlated with a high level of significance; > 5cr for 
log age and > 4cr for [Z/H] and [a/Fe]. Performing linear fits 
to the relations we find 

log age = (0.64 ± 0.12) log cr - (0.63 ± 0.26) , (2) 
[Z/H] = (0.40 ±0.08) log cr- (0.63 ±0.17) , (3) 
[a/Fe] = (0.20 ±0.06) log cr- (0.17 ±0.13) . (4) 

These fits are shown as red lines in Figure [S] While these 
results confirm the existence of the [Z/H]-cr and [a/Fe]-cr re- 
lations, the correlation of age with velocity dispersion found 
here is not so well accepted (see Section [T] for references). 

The behaviour of age and [ct /Fe] with velocity disper- 
sion is reminiscent of down-sizing l|Cowie et al.lll996l ) , where 
the typical mass of a star-forming galaxy increases with red- 
shift. Looking at the distribution of ages, we see that at all 
velocity dispersions there exist galaxies with very old stellar 
populations, but the age of the youngest galaxy at a given 
velocity dispersion increases with velocity dispersion. The 
peak epoch of star formation in more massive galaxies thus 
occurred at higher redshifts. 

Assuming that the scatter about the relations is Gaus- 
sian we find that the intrinsic scatter in log age is 0.20 ±0.02 
dex, in [Z/H] is 0.10 ± 0.02 dex and in [a/Fe] is 0.07 ± 0.01. 
The uncertainty in these scatters were determined via Monte 
Carlo simulations. Considering the mean observational er- 
rors are 0.14, 0.08, and 0.08 dex, respectively, the tight- 
ness of these relations is remarkable, and shows that the 
velocity dispersion is an excellent indicator of the SPPs 
in early-type galaxies, with more massive galaxies being 
older, more metal-rich and having shorter star-formation 
timescales than less massive galaxies. 

Table [5] lists the SPP-o" relations found in this work 
and in the literature. The literature estimates come from 
a sam ple of early-type g alaxies in high-density environ- 
ments l|Thomas et al|2005l). asarnple of red-sequence galax- 
ies in low-z clusters |Nelan et al.ll2005l ). a sample of early- 
type galaxies in high -density environments from the SDSS 
ISernardi et al. 20061). a sample of red-sequence galaxies in 
low-z clusters ( Smith et al.ll2006l), a sample of red- sequence 
galaxies drawn from the SDSS l|Graves et al.|[2007l ). a sam- 
ple of galaxies from the Shapley supercluster with cr > 100 



km s-i ((SmitheraD[2003), and a sample of red sequence 
galaxies in all environments. 

Overall, the slopes found here agree well with those 
found in the literature. The slope of our logage-cr relation 
is consiste nt with most e s timat es in the literature except 
for that of iThomas et al.] I 
steep as ours, and that of 



2005h. which is only a third as 



Bernardi etU] (|2006l ) who find 



a slope twice as steep. The existence of a correlation be- 
tween age and velocity dispersion is still being debated and 
this is reflected in the large variation in slopes found in the 
literature. A possible cause for this is the difficulty in accu- 
rately determining ages with the current set of models. The 
slope of our [Z/H]-a relation agre es well with those f ound 
in the literature, with exceptions of lGrave^ et all (|2007l ) and 
iThomas etakl 12010') who find a slope twice and 50 percent 
as steep, respectively. The slope of our [a/Fe]-cr relation is 
consistent with all of the literature values despite being the 
smallest. 



5.2 The relations in the cluster-outskirts sample 

The median values of the SPPs were found to be the same 
in both the cluster sample and the cluster-outskirts sam- 
ple, as were the SPP distributions. In order to make a true 
comparison, one that is free from any biases introduced by 
the two velocity dispersion distributions, we now compare 
the parameter-cr relations found in the cluster sample with 
those in the cluster-outskirts sample. 

Figure |5] shows the parameter-cr relations for the galax- 
ies in the cluster-outskirts sample (solid red lines) compared 
with those in the cluster sample (dashed black lines). Details 
of the fits are given in Table [S] 

None of the SPPs are found to be signiflcantly corre- 
lated with velocity dispersion in the cluster-outskirts sam- 
ple. Note that the parameter errors are larger in the cluster- 
outskirts sample than the cluster sample. This is due to the 
fact that the galaxies in the cluster outskirts were observed 
with 6dF, which has less wavelength coverage than 2dF, and 
therefore their SPPs were estimated, on average, from fewer 
indices. It is possible that these larger errors blur any cor- 
relation that may exist. 

It is also possible that the change in the number and 
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Table 5. Comparison of the SPP-o" relations from this study and the literature. 



Reference 


log a, 


'"e 


[Z/H] 


[a/Fe] 




Zero-point 


Slope 


Zero-point 


Slope 


Zero-point 


Slope 


This study 














Cluster 


-0.63±0.26 


0.64±0.12 


-0.63±0.17 


0.40±0.08 


-0.17±0.13 


0.20±0.06 


Cluster-outskirts 


-0.43±0.45 


0.54±0.20 


-0.24±0.40 


0.18±0.18 


-0.52±0.31 


0.33±0.14 


Thomas et al. (2005) 


0.46 


0.24 


-1.06 


0.55 


-0.42 


0.28 


Nelan et al. C2005~) 




0.59±0.13 




0.53±0.08 




0.31±0.06 


Bernardi et al. (2006) 


-1.72±0.31 


1.15 


-0.64±0.01 


0.38 


-0.54±0.01 


0.32 


Smith et al. f2006) 




0.72±0.14 




0.37±0.08 




0.35±0.07 


Graves et al. (2007') 




0.35±0.03 




0.79±0.05 




0.36±0.04 


Smith et al. f2007) 




0.64±0.12 




0.38±0.09 




0.36±0.07 


Thomas et al. ("2010) 


-0.11±0.05 


0.47±0.02 


-1.34±0.04 


0.65±0.02 


-0.55±0.02 


0.33±0.01 




combination of indices used to estimate the SPPs in the 
cluster-outskirts sample, when compared to the cluster sam- 
ple, is the reason why no correlations are found in the former. 
We therefore check this effect by re-estimating the SPPs of 
the cluster sample using only H/3, Mgb and Fe5335. Using 
these estimates we find that all three parameters are still cor- 
related with velocity dispersion and that the slopes of both 
relations are consistent within the errors. Therefore it is un- 
likely that the lack of correlations in the cluster-outskirts 
sample is a result of the change in the number of indices 
used. We would like to point out that a fit (using the un- 
restricted set of indices) was only accepted if H/3 was used 
and at least two other indices. See Paper I for more details. 

We show the results of comparing the two sets of es- 
timates in Figure [TOl For [Z/H] and [a/Fe] the agreement 
between the two sets of estimates is quite good, although 
there is an offset of 0.05 dex in [a/Fe] in the sense that the 
unrestricted set of indices gives larger values of [a/Fe]. There 
is also an offset of 1.3 Gyr between the two age estimates 
(but in the opposite sense to that of the [a/Fe] offset) and 
there is some scatter that increases with age. This is due 
to the fact that, in index space, the distance between lines 
of constant age decreases with increasing age so that small 
shifts in index strengths can lead to large changes in age 
estimates. According to a Spearman rank correlation test 
the two datasets are consistent. This implies that the two 
ways of measuring the stellar population parameters rank 



the objects in the same order and that there should be no 
effect on our results. 



6 THE SCATTER IN THE PARAMETER 
DISTRIBUTIONS IN CLUSTERS 

In part, the spread in the SPP distributions found for the 
cluster sample is due to a combination of our observational 
errors and the correlations with velocity dispersion. To de- 
termine the amount of intrinsic scatter in the distributions, 
over and above that caused by the observational errors, we 
ran a series of numerical simulations. The method used was 
as follows. 

Due to the need to create a large number of models 
containing a large number of galaxies, which need to have 
their SPPs converted to line-strengths (and vice versa), it 
is not feasible to estimate SPPs using all available indices. 
Therefore, for this exercise, we rely on the three indices used 
earlier in investigating the anti-correlation between age and 
[Z/H], i.e. H/3, Mg6 and Fe5335. 

With this combination of indices, galaxies in the clus- 
ter sample have an age distribution that is approximately 
exponential and [Z/H] and [a/Fe] have distributions that 
are approximately Gaussian (see Figure [5]). Therefore, we 
began by selecting a range of e-foldings for the exponential 
age distribution (r = 0.10, 0.15, 0.30, 0.55, 0.90, 1.60, 2.80, 




5.00, 8.50, 15.00 Gyr), and ranges of Gaussian scatters in 
[Z/H] (cr[z/H] = 0.03, 0.05, 0.08, 0.12, 0.20, 0.30, 0.50, 0.80, 
1.25, 2.00 dex), and [a/Fe] (cr[„/Fc] = 0.03, 0.04, 0.07, 0.10, 
0.15, 0.20, 0.30, 0.45, 0.70, 1.00 dex). Each range was cho- 
sen to be approximately evenly spaced logarithmically so as 
to have better resolution at small values. The medians of 
the cluster [Z/H] and [a/Fe] distributions were used as the 
means of their Gaussian distributions, and the exponential 
age distribution was truncated at 5 Gyr (because the models 
are less reliable at young ages) and at 14 Gyr (the age of 
the universe). 

A model was generated for each combination of r, 
o"[z/H], and o-^^/fc] consisting of 10,000 mock galaxies with 
ages, [Z/H], and [a/Fe] (parameter triples) drawn randomly 
from the specified distributions. If any of the mock galax- 
ies had a [Z/H] or [a/Fe] that fell outside the extent of the 
modelj^ then they were assigned the closest value from the 
models; this was not necessary for the ages since they were 
truncated at the values stated above. 

These pa rameter triples were converted, via the 
omas et al.l models, to the corresponding H/3, Mg6, and 
Fe5335 index values (index triples). This was achieved as 
follows. The models were divided up into cells defined by 
two adjacent values of each SPP, and for each mock galaxy 
the cell in which it fell was determined. Each cell is also de- 
fined by eight values of H/3, Mg6, and Fe5335, corresponding 
to their values at the eight vertices of the cubical cell. The 
index triple corresponding to the parameter triple is then 
calculated by taking the mean of each of these eight index 
values weighted by the distance in parameter space of the 
mock galaxy from each corner of the cell. 

These index triples were then perturbed by randomly 
drawing errors from an observed galaxy's index error distri- 
butions. In doing this we cycled through the observed galax- 
ies so that each observed galaxy had an approximately equal 
number of mock galaxies perturbed by its index error dis- 
tributions. The perturbed index triples were converted back 
to parameter triples. 

To ascertain how well the models fit the data a likeli- 



1 The lThomas et al.l ||2003|) models have 1 < t ^ 15 Gyr, -2.25 ^ 
[Z/H] sC 0.67 dex, and 0.0 [a/Fe] ^ 0.5 dex. 



hood statistic was used. The SPP space was divided up into 
bins and the probability of each bin containing a galaxy 
was calculated. The age bins were 0.5 Gyr in width and 
ranged from and 15 Gyr, the [Z/H] bins were 0.125 dex 
in width and ranged from —2.2875 to 0.7125 dex, and the 
[a/Fe] bins were 0.025 dex in width and ranged from 0.0 to 
0.5 dex. This resulted in 30 age bins, 24 [Z/H] bins, and 20 
[a/Fe] bins. Before the likelihood statistic was calculated, 
the model was lightly smoothed to minimise the number of 
bins with zero probabilities; a 3D Gaussian, with a = 2/3 
bins for all dimensions, was used as the smoothing kernel, 
and bins within 3a, in each dimension, were used in calculat- 
ing the value of the bin probability in the smoothed model 
(i.e. two bins either side in each dimension were used). The 
likelihood statistic is obtained by determining which bins 
the observed galaxies lie in and summing the logarithm of 
the bin probabilities. 

Monte Carlo simulations were used to estimate the 
probability of obtaining a likelihood statistic larger than the 
observed one by chance. There are 103 galaxies with mea- 
surements of the indices required to estimate the SPPs in the 
cluster sample, so we simulate 103 galaxies in an identical 
manner to the model and calculate the maximum likelihood 
statistic. This is done 10,000 times. The fraction of simula- 
tions with likelihood statistics larger than the observed value 
are recorded and 3D confidence contours are generated. 

These 3D confidence contours are shown in Figure [TTl 
which presents the likelihood of our dataset being repre- 
sented by a model with the specified e-folding of the expo- 
nential distribution of ages and specified RMS scatters in 
the Gaussian distributions of [Z/H] and [a/Fe]. While Fig- 
ure [TT] provides an overall picture of the acceptable scatters 
in the SPPs, more detailed information can be revealed by 
taking slices through the 3D confidence contours. Figure [T2l 
shows slices taken parallel to the age- [Z/H] plane at val- 
ues of [a/Fe] corresponding to those used in generating the 
models. 

These simulations provide us with a great deal of infor- 
mation regarding the intrinsic scatter in each of the SPPs. 
Firstly, it is difficult to constrain the e-folding of the expo- 
nential age distribution. This is due to small uncertainties 
in line-strengths translating to large changes in the age es- 
timate combined with relatively large observational errors 



16 Craig Harrison et al. 




Figure 11. Two views of the 3D confidence contours of the e-folding of the exponential distribution of ages (r) and the RMS scatters 
in the Gaussian distributions of [Z/H] and [«/Fe] (o-[2/h] s-i^d ^[a/Fe])i generated from the numerical simulations described in the text. 
The blue, green, and red surfaces are the Icr, 2a, and 3cr confidence contours. 
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in the H/3 index. Although models with small values of e- 
folding are not ruled out, we find that our data is most con- 
sistent with the model having r = 900 Myr. Secondly, small 
scatters in [Z/H] (o"[z/h] < O.f dex) are strongly ruled out, 
as are large scatters (o"[z/h] > 2.0 dex). We find the model 
with cr[z/H] = 0.3 dex is the most consistent with our data. 
Finally, large scatters in [a/Fe] (fj^/pcj > 0.3 dex) are also 
strongly ruled out. Models with low scatters in [a/Fe] are 
not ruled out, but we find that the most consistent model 
is that with a^a/pc] = 0.07 dex. The model that gives the 
maximum likelihood (i.e. r ~ 900 Myr, cr[z/H] ~ 0.3 dex, 
and cr^a/Fc] = 0.07 dex) is marked by a cross in Figure [121 

Since the SPPs were found to be correlated with velocity 
dispersion, and since the galaxies in the cluster sample have 
a range of velocity dispersions, it is understandable that we 
should detect an intrinsic scatter in the SPP distributions. 
The degree to which the intrinsic scatter in the parameter 
distributions is attributable to trends with velocity disper- 
sion can be estimated by comparing the scatter expected 
given the velocity dispersion distribution in the cluster and 
the intrinsic scatter about the parameter-a relations. 

We use the parameter-cr relations to convert the galax- 
ies' velocity dispersions into parameter values. The scat- 
ter in these values is then calculated; the intrinsic scatter 
about the relation (as found in Section 15.11) is then added 
in quadrature to the scatter due to the relation with veloc- 
ity dispersion, and the result is compared to that obtained 
from the simulations described above. We do this only for 
[Z/H] and [a/Fe], which were found to have approximately 
Gaussian distributions. 

The scatter expected in [Z/H] on this basis is 0.21 dex 
which is comparable to the intrinsic scatter of cr[z/H] ~ 0.3 
dex estimated above. The expected scatter in [a/Fe] is 0.06 
dex which is very close to the estimated intrinsic scatter 
c"[a/Fc] ~ 0.07 dex. It appears then that the intrinsic scatter 
in both the [Z/H] and [a/Fe] distributions can almost en- 
tirely be accounted for by the parameter-a relation and the 
intrinsic scatter about it. 

For [Z/H], the scatter due to the correlation is 0.16 dex 
while the intrinsic scatter about the relation is 0.13 dex; 
thus, most of the intrinsic scatter in the [Z/H] distribu- 
tion comes from the correlation with velocity dispersion. For 
[a/Fe], the scatter due to the correlation is 0.03 dex while 
the intrinsic scatter about the relation is 0.05 dex; thus most 
of the intrinsic scatter in the [a/Fe] distribution comes from 
the intrinsic scatter about the [a/Fe]-cr relation. 



7 DISCUSSION 

The process by which galaxies form and the factors with the 
greatest influence on their evolution are unresolved issues. 
In models of galaxy forma tion based on hierarch ical merging 
in a ACDM universe (e.g. iDe Lucia et al.|[2006l ) the merger 
history of an early-type galaxy can fall anywhere between 
two extremes: mergers could occur early-on and rapidly with 
subsequent passive evolution — th e revised monolithic col- 
lapse scenario (e.g. iMerlin fc Chiosi 2006 ); or the galaxy 
could experience a more prolo nged history o f mergers — the 
extended merging scenario (e.g. lToomrj|l977l ). The environ- 
ment of a galaxy and its mass are the two most influential 
factors on galaxy evolution but their effects need to be to 



disentangled to gain an insight into the relative importance 
of each. 

For cluster galaxies, the numerical simulations per- 
formed here show that the intrinsic scatters in [Z/H] and 
[a/Fe], and the e- folding of the exponential distribution 
of ages, are all quite small, as would be expected from 
a rapid formation at high redshift followed by passive 
evolution. Additionally, positive correlations are found 
between all three of the SPPs and velocity dispersion. 
The correlations fo und for both [Z/H] and [a/Fe] confirm 
well- kn own tre nds (Trager et al.' '2000a'; 'Kuntschner et akl 
20011: [Proctor & Sansom 2002; Thomas ct al. 200a 
Caldwell et all 2003; Me hlert et al.li2003l: iNelan et al.ll2005l: 
Thomas et'alT i2005. : Bernardi et al.l l2006l: ICallazzi et all 



20061: iKelson et al 



Thomas et akl l2010h 



l2006l : ISanchez-Blazauez et al.l l200a 

However the correlation found 



for age is not widely recognised; some authors find 

to be positively corre- 



lated JProctor & SansomI 2002 


: Proctor et al.l 2004a 3: 


Nelan et al. 12005: Thomas et al.l 


20051: Gallazzi et al. 2009: 


Bernardi et all l2006l: iThomas et al.l 


20101) while others 


find they are not correlated at all 


Trager et all l2000al: 


Kuntschner et al.ll200ll: iTerlevich & Forbes! 20021). Not onlv 



do we find that age is positively correlated, but also that it 
is the most significantly correlated out of the three SPPs, 
being significant at the 5a level. These correlations are 
qualitati vely consistent with more recent semi-analytic 
models (|De Lucia et aL 20061'). These m odels incorporate 
feedback from AGN (|Croton et all 120061 ') , which heats the 
available gas, preventing further star formation in massive 
galaxies and bringing the predicted scaling relations into 
agreement with those observed. 

There is a level of homogeneity to early-type cluster 
galaxies evidenced in the similarity of their line-strength 
and SPP distributions and the consistency of the correla- 
tions between these quantities and velocity dispersion that 
is found between the clusters studied here. That this homo- 
geneity is maintained amongst clusters of varying richness 
and morphology indicates that differences in the cluster en- 
vironment have relatively little effect on the stellar popula- 
tions of early-type galaxies. The dominant factor appears to 
be mass. 

While mass plays a large role in determining the SPPs 
in early-type galaxies, it does not do so completely. The size 
of the role played varies; while [Z/H] is almost entirely a 
function of velocity dispersion, [a/Fe] is less so. The intrin- 
sic scatters in [Z/H] and [a/Fe] in the cluster galaxies were 
found to be almost entirely accounted for by the scatter pro- 
duced by the correlations with velocity dispersion and the 
intrinsic scatter about these relations. 

Galaxies in low-density environments are observed 



to be, on averagej_ 



2 Gy r youn g er than galajc- 



ies i n clusters |Trager et alT 2000al: Poggianti et al.l 



20011: iKuntschner et alj |2002|: iT crlevich fc Forbes! 1200"^ 
Caldwell et"al] !2003!: !Proctor et aL i2004bi: !Denicol6 et al] 



20051: 



Thomas et al 



Bernardi et al.l 



2006! : 



Sa nchez-Blazguez et al.! !2006! ). The situation is not so 
clear with regards to differences in [Z/H] and [a/Fe]. 
Recent models predict that galaxies in denser environ- 
ments are older and m ore metal-rich than isolated galaxies 
(|De Lucia et al.l!2006! ). We do not find significant differences 
in the line-strength distributions or the SPP distributions 
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between the cluster-outskirts sample and the cluster sam- 
ple. However, the size of our SPP errors make it difficult 
to detect the small differences in age reported by other 
authors and any difference may be masked by the f .3 Gyr 
offset to older ages found for galaxies that their parameters 
estimated from fewer indices; such galaxies are found 
predominantly in the cluster-outskirts sample. We do find 
that the tight correlations with velocity dispersion of both 
the line-strengths and SPPs that are found in clusters are 
weaker in the cluster outskirts, suggesting that the modes 
of formation in the cluster outskirts are more varied than 
those in the cluster cores. 

Similar results were found bv lThomas et al.1 (|2010l ) who 
conclude that the formation and evolution of early-type 
galaxies are relatively insensitive to environment and are in- 
stead driven by self-regulation processes and their intrinsic 
properties such as mass. However, our result that the SPP 
trends with velocity dispersion are weaker in the cluster- 
outskirts tlmnjii_th«£lu with these results 
(see also iBernardi et"aL I ll998l . I2OO6D . The dip in ages found 
near the cluster virial radius (see Figure [Tjl suggests that 
movement into the cluster environment may induce a burst 
of star formation and, combined with the down-sizing evi- 
dent in Figure [S] that this burst may occur preferentially 
in lower m ass galaxies (si milar to the rejuvenated popu- 
lation from iThomas et al.l ). Presumably, our two samples 
(cluster and cluster-outskirts) are combined in their single 
high-density sample, which may result in the masking of this 
trend. Our r e sults are broadly consistent with the results of 
ISmith et al] (I2OO6I') who find a variation of the SPPs with 
cluster radius. Thomas et al.l find that environment becomes 
more of an irifiuenc e moving to less massive galaxies (see also 
iHaines et aLll2006l ). 

Our results can be best explained by the revised mono- 
lithic collapse model, but is also consistent with the ex- 
ten ded merging model if the mergers are dissipationless 
fsee lDe Lucia et aLllioO^ L lPipino et al.l (|2009l ) have a model 
of elliptical galaxy formation that implements a detailed 
treatment of chemical evolution, and they find that the 
correlation between [a/Fe] an d velocity dis persion is much 
less steep than that observed. iPipino et al.l find that AGN 
quenching can help to improve the agreement, but this 
worsens their mass-metallicity relation compared to obser- 
vations. They suggest that both relations can be repro- 
duced provided the formation of all spheroids happens quasi- 
monolithically, i.e. the formation of the stars occurs in sub- 
units and that this star formation and the assembly of the 
sub-units into an elliptical galaxy occur at the same time 
and in the same place. 

The star-formation rates in galaxies are found to in- 
crease with increasing distance from the cluster centre, and 
to converge to the mean field rate at distances greater than 
~ 3J?vir time s the virial radius (Lewis ct al. 2002, see also 
iGomez et al.l (1200311 ). We find evidence that there is a dip 
in the mean ages of galaxies just inside the cluster virial ra- 
dius, possibly due to secondary star formation that reduces 
the luminosity-weighted mean age. This could be evidence 
that in-falling galaxies, upon reaching the virial radius of a 
cluster, undergo a burst of star formation triggered by the 
dense intra-cluster medium. If such a burst only amounted 
to a small fraction of the galaxy's total mass then its effect 
on the integrated light would be short-lived and the galaxy 



would rapi dly return to appearing old and red. Similar to the 
findings of iLewis et al.l . we find that this decrease in age ac- 
tually begins at > 3i?vir, suggesting that the influence of the 
cluster e nvironment extends to large distances. Whereas the 
result of lLewis et al.l was based mainly on late-type galaxies, 
the remarkable thing about this result is that it applies to 
early- type galaxies. 

There are two sources of bias that we must check have 
not influenced the above results: the number of indices used 
to derive the SPPs and and the aperture corrections. It is 
possible that the lack of correlations with velocity disper- 
sion for the SPPs in the cluster-outskirts sample is due to 
the fact that, on average, fewer indices were used to esti- 
mate them than those in the cluster sample. This results in 
larger SPP errors and might potentially explain the lack of 
correlations with velocity dispersion. However, as we show 
in Section [5. 2 1 using only three indices to estimate the SPPs 
does not signiflcantly alter their distributions compared to 
those derived from an unrestricted set of indices. 

It is also possible that mass segregation within the clus- 
ter could cause the weakening of the line-strengths in the 
cluster outskirts, if that is where less massive galaxies are 
preferentially found. This is a result of the same aperture 
correction being applied to a galaxy irrespective of how mas- 
sive it is and despite the fact that the line-strengths be- 
ing corrected were measured within differing effective radii. 
We find that our two samples (cluster and cluster-outskirts) 
have a similar shape to their magnitude distributions but 
that the galaxies in the cluster-outskirts are on average 
~ 1 mag brighter, due to the facts that there are more galax- 
ies in the cluster sample and that we target the brightest 
galaxies. So the correlation of weak line-strengths for less 
massive galaxies cannot be the reason for the weakening of 
the trends with velocity dispersion in the cluster-outskirts 
galaxies. We see no reason, therefore, not to believe that the 
weakening of the relations in the cluster-outskirts is real. 

The overall picture that emerges from this study is 
as follows. Early-type galaxies in clusters form a homoge- 
neous class of objects that form in a process similar to a 
revised monolithic collapse and whose stellar populations 
are largely determined by their velocity dispersion (mass) 
and are relatively unaffected (at least differentially) by the 
cluster environment. The more massive they are the older 
and more metal-rich they are, and the shorter their star- 
formation timescales. The stellar populations of early-type 
galaxies in the outskirts of clusters, in contrast, appear less 
influenced by mass because, due to the varying environments 
they formed and evolved in, their evolutionary histories are 
more varied and this causes the correlations with mass to be 
less significant. A galaxy, especially a less massive one, that 
is falling into a cluster will, upon nearing the virial radius, 
undergo a burst of star formation. Once these stars have 
faded and ceased dominating the integrated light the galaxy 
will appear indistinguishable from those in the cluster core. 



8 SUMMARY 

In summary, we have measured velocity dispersions, red- 
shifts and absorption-line strengths for a magnitude-limited 
(bj ^ 19.45) sample of galaxies drawn from four clus- 
ters (Coma, A1139, A3558, and A930 at (2} = 0.04) and 
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their surrounds (extending to ~ lOiivir)- Using the fully- 
cahbrated absorption- hne indices co u pled with the stellar 
population models of iThomas" et al.l l|2003l ). we have esti- 
mated ages, [Z/H], and [a/Fe] for 219 galaxies. We have 
used these data to investigate the effects of mass and en- 
vironment on the stellar populations of early-type galaxies 
and our results can be summarised as follows. 

(i) For galaxies in the cluster sample, all indices are pos- 
itively correlated with velocity dispersion, with the excep- 
tions of H/3 and H[/3g, which are negatively correlated, and 
Fe5406, which is uncorrelated. 

(ii) Only Mg6 and H/3g are correlated with velocity dis- 
persion in the cluster-outskirts sample. The slopes of these 
two relations are consistent with those found for the cluster 
sample. 

(iii) The cluster cores are relatively free from young galax- 
ies and from galaxies that have experienced recent star for- 
mation. These galaxies are more commonly found outside 

-RAbcll- 

(iv) The stellar populations in clusters form a homoge- 
neous population. Despite the fact that our sample was 
drawn from four clusters spanning the ranges of Bautz- 
Morgan classifications and Abell richness classes, the line- 
strength-a relations, parameter-cr relations, and SPP distri- 
butions are consistent between clusters. 

(v) There is no difference between the line-strength dis- 
tributions and SPP distributions in the clusters and their 
outskirts. 

(vi) The SPPs in clusters are correlated with velocity dis- 
persion, suggesting that more massive galaxies are older, 
have shorter star-formation timescales, and are more metal- 
rich. These correlations are found to be weaker in the cluster 
outskirts. 

(vii) For galaxies in the cluster sample, the e-folding time 
of the age distribution is 900 Myr, the intrinsic scatter in the 
[Z/H] distribution is 0.3 dex, and the intrinsic scatter in the 
[a/Fe] distribution is 0.07 dex. These latter two intrinsic 
scatters can almost entirely be accounted for by the corre- 
lations with velocity dispersion and the scatter about the 
relations. We conclude, therefore, that the mass of a galaxy 
plays a major role in determining its stellar populations. 

Further high quality observations of galaxies at higher 
redshifts will allow the development of a consistent model 
of early-type galaxy formation at all masses and in all 
environments. 

The authors thank the anonymous referee for a careful 
reading of the paper and for the many helpful suggestions 
that improved it considerably. 
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